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PREFACE 


This report was prepared under contract NAS 1-12579, "Expansion 
and Extension of the SBOOM Computer Program. " The study was 
carried out during the period from July, 1973 through December, 
1973. The study was funded tiy the National Aeronautics and Space 
Administration, Langley Research Center, Space Systems Division. 

The study effort resulted in the development of two new computer 
programs, one for generating and maintaining a data base of near- 
field pressure signatures and the other for predicting sonic boom 
as a result of overflight of shuttle type reentry vehicles. The 
study results extend the work perfoimned under contract NAS 2-6147 
to NASA Ames Research Center in which the basic method of predict- 
ing and optimizing shuttle trajectories based on sonic boom con- 
straints was developed. Both contracts employed a pressure 
signature extrapolation technique and wind tunnel measurements 
developed by Ames Research Center. 

The data base management system developed for the original contract 
and used extensively for two dimensional sonic boom prediction 
proved inadequate for the three-dimensional requirements of the 
present contract. A new data management system was developed which 
is versatile enough to handle present and future needs of the 
sonic boom methods employed, ' 

The report is presented in two volumes: . 

Volume I Method and' Results 

Volume II Data Base Construction 

The first volume describes the method employed for estimating 
ground overpressures from wind tunnel measurement. Some results 
are presented and the use of the computer program is described. 
Volijme II describes the data management system employed in the 
data base construction and maintenance. A separate computer pro- 
gram developed for this purpose is also described. 
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PREDICTION OF SONIC BOOM FROM EXPERIMENTAL NEAR-FIELD 


OVERPRESSURE DATA. VOLUME I - METHOD AND RESULTS. 

BY C. R. Glatt/ D. S. Hague and S. J. Reiners 
Aerophysics Research Corporation 

SUMMARY 


A technique for prediction of ground level sonic boom overpressures 
generated by the flight of space shuttle vehicles is presented. 

The method is in the form of a digital computer program which con- 
tains the following key elements: 

1. A data base of typical near-field space shuttle orbiter 
pressure signatures generated by extensive wind tunnel 
testing at specified flight conditions. 

2. An interpolation procedure which operates on the near- 
field space shuttle signature data base to produce the 
near-field space shuttle pressure signature for any flight 
condition. 

3. A trajectory analysis module capable of generating the 
space shuttle flight path. 

4. A method for estimating ground level pressure signatures 
corresponding to each flight condition encountered along 
the space shuttle trajectory. This method is based on 
extrapolation of a near-field pressure signature to sea 
level conditions. 

5. An optimization procedure by which the sea level over- 
pressures generated are reduced by trajectory perturbation. 

6. A method of automatically generating graphical representa- 
tions of measured and estimated pressure signatures. 

The technique is discussed in detail and examples are presented 
which illustrate flight paths which generate sea level overpressures. 
The experimental pressure signatures used to construct the near- 
field signature data base and the procedure for extrapolating a 
near-field signature to sea level were developed by NASA personnel 
prior to this study. The near-field signature data base and inter- 
polation procedure were developed as the result of prior government 
sponsored research. The present report describes the extension 
of the method to encompass roll angle predictions as well as Mach 
number and angle of attack. 



INTRODUCTION 


A sonic boom computer program (SBOOM) based on the prediction 
technique of Thomas (reference 1) is presented. The prediction 
technique calculates far-field overpressure from near-field 
pressure signatures measured in the wind tunnel. The program is 
used as an independent program for the determination of sonic 
boom characteristics of a specified vehicle and environment for 
input flight conditions or the program can be linked to a trajec- 
tory prediction program through the use of the linking system 
described in Reference 2. 

Wind tunnel data for a space shuttle delta wing orbiter configura- 
tion has been digitized and stored as a pressure signature data 
base using a newly developed access and retrieval system which 
provides unlimited expansion as more signature data is acquired. 

A program, GETTAB, has been written for the purpose of storing 
measured pressure signatures in the data base. The program also 
performs certain mapping of known signatures into signatures for 
flight regions where data is unavailable. The data base is accessed 
by the SBOOM program and pressure signatures for input flight 
conditions can be estimated using the geometric similarity rule 
developed. The data base and procedure developed estimates pressure 
signatures in the flight regime. 

Angle of Attack (10 < a < 45) 

Mach Number (1.2< M < 10) 

Roll Angle (0 < <j) < 180) 

The flight regimes span the normal range expected during reentry 
but does not cover all conditions ' expected during launch. Tests 
at Ames Research Center indicate that plume effects invalidate the 
data under launch conditions unless a "plume factor" is employed. 

A research effort to generate launch data including plume effects 
is now underway at Ames. 

The original SBOOM program was incorporated into the Trajectory 
Optimization Program (ATOP) of references 3 through 5 as an 
independent link to form the Sonic Boom Trajectory Optimization 
Program. In operation, the ATOP program generates a special file 
of data containing the trajectory information needed by the SBOOM 
subprogram. SBOOM then interrogates the special file of informa- 
tion in a repetitious manner generating sonic boom overpressure 
data at each point along the trajectory. The SBOOM program reported 
here has been isolated as an independent program. As such, SBOOM 
can obtain trajectory information from any trajectory program and 
still access the data base of pressure signatures. Options are 
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available for reading trajectory tapes from ATOP and the Program 
to Optimize Simulated Trajectories (POST, unpublished NASA pro- 
gram) . 

The combined SBOOM/POST/AESOP simulation illustrated in figure 1 
can be used in optimizing both pitch-plane and roll-pitch-plane 
trajectories for the shuttle vehicle with sonic boom constraints. 
The AESOP parameter optimization program of reference 6 is 
employed for this purpose. The ODIN Executive System of refer- 
ence 2 is used to link the independent programs POST, SBOOM and 
AESOP to optimize trajectories based on sonic boom constraints. 

An additional feature of the SBOOM and GETTAB programs is the 
ability of these programs to generate plots of the pressure 
signature data from the data base. This is done by an interface 
with a separate plotting program. 

This report is presented in two volxames. Volume I contains a 
discussion of the SBOOM program and results obtained from the 
use of the data base to estimate near-field pressure signatures. 
Voliime II presents a discussion of the GETTAB program, the method 
of access and retrieval, and plots of the entire pressure signa- 
ture data base. 
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FIGURE 1 ILLUSTRATION OF THE USE OF SBOOM WITHIN ODIN 
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DISCUSSION OF THE PROGRAM 
Sonic Boom Prediction 


The computer program described in this report uses the 
method of Thomas (reference 1) for extrapolating near- 
field signatures out to the far field, without use of 
a F-function to account for nonlinear waveform distortion. 
Effects of aircraft acceleration and atmospheric tempera- 
ture, pressure and wind gradients are included in the 
theory. The approach used is to describe the waveform 
of the sonic boom wave by several waveform parameters and 
then to obtain equations for the parameters as functions 
of time. This approach has the advantages that (1) the 
theory is simpler and more intuitive than the Whitham 
theory, (2) it provides a more convenient method for 
extrapolating experimental signatures because the signa- 
ture is dealt with directly, rather than through the use 
of a F-function and (3) shock locations are determined 
by a much neater method than the classical area balancing 
technique used in F-function extrapolations. 

To describe the waveform of the sonic boom wave at any 
instant of time, we approximate the waveform by an arbi- 
trary number of linear segments and define the waveform 
parameters Ap., m. and X. of each segment as follows: 

Ap^ is the pressure rise across the shock at the juncture 

of segments i and i - 1. Often there will be no shock at 
the juncture, in which case Ap^^ is zero, m^ is the slope 

of segment i, which may be positive or negative. Finally, 
is the length of segment i. A completely general wave- 
form can be described using these waveform parameters. 

To determine the waveform parameters as functions of time, 
it is assumed that the time rate of change of any waveform 
parameter can be obtained by superposition of the rate of 
change assuming the wave propagates as a linear, nonplane 
wave and the rate of change assuming the wave propagates 
as a nonlinear, plane wave; it is, for example: 



nonlinear 

Tionsdano 




linear 
nvni'Ianc 


d^Pi \ 
dt ) 


nonlinear 
plane ' 
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The linear wave term accounts for the affects of changing 
ray tube area and changing atmospheric properties. The 
nonlinear, plane wave term, which is also influenced by 
the atmospheric properties, accounts for the nonlinear 
distortion of the waveform. Details of the waveform 
parameter method are given in reference 1. 


Pressure Signature Data Base 

The procedure for determining sonic boom overpressures 
on the ground produced by vehicles flying at supersonic 
speeds is to define the near-field pressure signatures 
and then to extrapolate these signatures to the far-field 
(ground) . Experience has shown that the best estimates 
of ground overpressure can be obtained by resorting to 
experiment rather than theory to determine near-field 
pressure signatures. 


The sonic boom pressure signature data base is a collection of 
the results of several wind tunnel experiments in which near- 
field pressure signatures are measured for a range of anticipated 
flight conditions. These experiments have been conducted for a 
delta-wing shuttle vehicle. The Mach niombers ranged from 1.2 to 
10.2. The angles of attack and roll angles ranged from 
0 to 60 degrees and 0 to 180 degrees respectively. Even 
though current booster and orbiter configurations may 
differ from the one being investigated, the pressure 
signatures for the shape being tested are reasonably 
representative of the ones of interest provided the angle 
of attack is not small. This is borne out of comparison 
of results of straight and delta wing shapes discussed 
in reference 7. 

Models having the shape shown in figure 2 have been tested 
in the Ames 3.35 by 3.35 meter (11 by 11 foot) 2.74 by 
2.13 meter (9x7 foot) and 2.44 by 2.13 meter (8x7 foot) 
wind tunnels, and the jet propulsion laboratory .508 meter 
(20 inch) supersonic and the 0.533 meter (21 inch) hyper- 
sonic wind tunnels. The various angles of attack and roll 
angles were obtained by rotating the model and sting 
assembly relative to the pressure measuring equipment. 

The models were mounted on a linear actuator which permitted 
them to be translated longitudinally in the wind tunnels 
relative to the fixed pressure measuring equipment. Flow 
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field pressures were detected by two degree included 
angle conical static pressure probe mounted on the wind 
tunnel wall which was connected to a capacitance type 
pressure transducer. Measured near-field sonic boom 
overpressures are presented in reference 7. A 
summary of the test conditions which resulted are pre- 
sented in figure 3. 


The wind tunnel test results are converted to digital form 
so that the pressure signatures can be used directly in the 
SBOOM computer program. To insure the validity of the data 
transcribed from graph to digital data, CALCOMP plots of 
the digitalized data were generated using an independent 
plotting program. These plots were drawn to the exact 
scale of those presented in reference 7. They 
were carefully compared with the original plots by over- 
laying one on the other. 

After the comparison of the digitalized data with the 
wind tunnel results, the data was transferred directly 
into the pressure signature data base using the auxiliary 
storage and retrieval program GETTAB described in Volume II. 
These data are available on data cell and may be used by 
simply attaching the proper file when executing the SBOOM 
program. 


Interpolation Procedures 


The use of stratified data represented by the wind tunnel 
test conditions required the development of some inter- 
polation procedures based on the geometric similarity of 
the pressure signatures. 


Interpol a ti on. - The program contains an interpolation 
rule based on geometric similarity of the actual pressure 
signatures measured in the wind tunnel. This approach, 
illustrated in figure 4, considers each near-field pressure 
signature to be comprised of four segments as follows: 
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FIGURE 4 ILLUSTRATION OF THE GEOMETRIC 
SIMILARITY RULE. 
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First significant overpressure to maximum overpressure. 

2. Maximiom overpressure to first following zero overpressure. 

3. First zero overpressure (after max DP/P) to minimum over- 
pressure. 

4. Minimum overpressure to last significant point in signature. 

NOTE: A Trailing Zero Pressure Point Should be Added to all 

Signatures where They are Omitted. 

Corresponding signature segments are mapped into each other using 
linear interpolation within the segments and a variable which spans 
the range 0 to 1 within each segment. Some possible signature 
types are illustrated in figure 5 which illustrates the diversity 
of signature type within the data base. The method developed will 
uniquely map any signature of figure 5 into any other signature 
of figure 5. Some interpolations at the mid point between pairs 
of idealized signatures are illustrated in figure 6. Some con- 
tinuous mappings of one signature into another are presented in 
figure 7 using signatures actually contained within the data base. 
Figure 7 indicates which pair of the signatures contained in ref- 
erence 7 are used and whether an interpolation or extrapolation 
is used. 

Interpolation between given signatures proceed linearly between like 
segments. An interpolation for the corresponding h/Z value (distance 
from the model to the point of overpressure measurement normalized 
to the model length) is also determined simultaneously with the 
signature interpolation by linear interpolation. Signatures produced 
in this manner have characteristics geometrically similar to the given 
signatures and the method provides a fairly accurate interpolation 
procedure. 


Comprehensive testing of the technique included interpolation 
between known signatures for a known signature. These tests showed 
good correlation between the known and the interpolated signatures 
as illustrated in figure 8. Typical results of extrapolating an 
interpolated near-field signature to the ground and an original 
data base signature to the ground are presented in figure 9. 
Interpolation by the geometric similarity rule always interpolates 
first in the roll angle plane, then the Mach plane and finally 
the angle of attack plane. This produces better results than any 
other order of interpolation. Figure 10 illustrates the inter- 
polation procedure and figure 11 shows some typical results. 

Near-Field Extrapolation 

Since data is not always available for interpolation, a polynomial 
extrapolation technique was developed for generating pressure 
signatures outside the range of the test data. The technique 
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EXTRAPOLATION OF INTERPOLATED NEAR-FIELD SIGNATURES 
TO SEA LEVEL (MACH NUMBER INTERPOLATION). 
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DATA BASE PRESSURE SIGNATUgES 
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FIGURE 10 ILLUSTRATION OF PRESSURE SIGNATURE 
INTERPOLATION PROCEDURE. 
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ALPHA (A) = 10 DEGREES. 
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FIGURE UB BANK ANGLE (P) INTERPOLATIONS AT MACH (M) = 3.98 AND 5.96 

ALPHA (A) = 10 DEGREES. 
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FIGURE IID ANGLE OF ATTACK (A) INTERPOLATION. 


developed used the input stations {h\/l) of the given curve and 
scales the AP/P values of the given curve by the following factor: 

(Ml/M2)^ 

where Ml is the Mach n\imber of the given curve, M2 is the Mach 
number of the desired curve, and k is a constant equal to 2 for 
extrapolating high Mach number signatures to higher Mach numbers 
and 1.5 for extrapolating low Mach number signatures to lower 
Mach numbers. No methods have been developed for extrapolating 
in the angle of attack or roll angle planes. 

Tests of the extrapolation technique included extrapolation of 
given experimental data base signatures for other data base experi- 
mental curves. Results of these curves showed the extrapolation 
technique used was quite accurate in the Mach plane, figures 12 
and 13. 


Enriching the Near-Field Signature Data Base 


The near-field interpolation and extrapolation procedures developed 
were used to produce an enriched near-field data base of signatures. 
In this way the amount of uncontrolled extrapolation of near-field 
experimental signatures during a trajectory analysis is minimized. 
Figure 14 illustrates the flight conditions for which wind tunnel 
measurements of pressure signatures were available. Low and high 
angle of attack data was very sparse so were omitted from the final 
data base. Mach 2.7 data was omitted because the interpolation 
method for augmenting this data produces the same results as those 
obtained from the SBOOM program. The flight conditions stored in 
the final data base are shown in figure 15. Estimations of pressure 
signatures outside this region will generate the boundary signature. 
The arrows in figure 15 indicate the source of the data for the 
pressure signature augmentations. 

Maximum ground overpressures obtained by a far-field extrapolation 
of the enriched data base of near-field signatures is summarized 
in figure 16. The data ranges spanned in the original tests of 
reference 7 are indicated by solid symbols. The smooth 
variation of maximum ground overpressure with angle of attack and 
Mach number appears to verify the procedure. All results in 
figure 16 are obtained by a far-field extrapolation of the stored 
near-field signatures from unaccelerated flight at 30480 meters. 
Near-field signatures are stored at each point indicated by symbols 
in the enriched data base. 
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FIGURE 12 
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FIGURE 16 RESULTS OBTAINED FROM MODIFIED PROCEDURE. 


RESULTS 


The technique for estimation of space shuttle ground overpressures 
described in the preceding section is summarized in figure 17. 
Measured near-field wing tunnel signatures are collected into a 
data base of known signatures as a function of: 

Mach nijmber, M 

Angle of attack a 

Roll angle, (p 

Given any two signatures in the data base, an interpolation pro- 
cedure based on geometric similarity produces an intermediate 
signature. By successive interpolation the signature correspond- 
ing to any point P in the (M, a, 4>) space can be found from the 
signatures stored at the discreet points “ j / The inter- 

polated signature, which corresponds to a flight condition generated 
by the trajectory is then extrapolated to the ground using the 
Thomas waveform parameter method. This procedure is repeated for 
all or a selected set of points along the trajectory and a measure 
of the ground overpressure is calculated for the trajectory. This 
measure may be any of several functions including: 



where the r points correspond to the selected set of trajectory 
points R^ in number. 



Where t^ and t 2 define the trajectory period of interest, t^^< t < t 2 , 

and ^p is the acceptable overpressure. Successive trajectories are 
then generated by systematically perturbing the vehicle control 
variable history, usually angle of attack or pitch angle. The 
systematic control variable perturbations are generated by the 
application of multivariable search technique, reference 6, which 
generally fall into one of several classes of search, for example: 

1. One parameter at a time. 

2. Systematic multi-parameter, such as steepest-descent or 
second order. 

3. Randomized. 
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SCHEMATIC 


FIGURE 17 SCHEMATIC OF SONIC BOOH ESTIMATION TECHNIQUE. 
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4. Acceleration, which improves the results previously- 
obtained by searches of Types 1, 2 or 3. 

The objective of these searches is to define the control history 
which minimizes or F 2 along the trajectory and thus the degree 

of objection to overland flight by space shuttle vehicles. It 
should be noted that, strictly speaking, the flight path optimiza- 
tion problem posed is a complicated problem in the variational 
calculus. The manner in which such a problem can be transformed 
into simpler multivariable search class by control parameteriza- 
tion has been discussed in detail in reference 6. Basically, the 
continuous control history, a(t), is replaced by a parameteriza- 
tion, a^(T^). Between any two time points, T^ and T^ + 1, an 

interpolation rule is used to define the instantaneous control. 

The finite set of parameters, a^(Tj^), then form the basis for 

applying multivariable search techniques to the trajectory optimiza- 
tion problem. 

Behavior of the ground oveiqsressure measure F^ along a typical 

shuttle ascent path is presented in figure 18. Here, the trajec- 
tory is typified by the Mach altitude schedule flown with the F^^ 

altitude schedule superimposed. The path illustrated is a high 
dynamic pressure path. Ground overpressure peaks early in the 
flight with the vehicle in a low altitude/low supersonic Mach num- 
ber condition. For the vehicle and path combination chosen, signi- 
ficant overpressures approximating 1 psf are still propagated to 
the ground as the vehicle passes through 100,000 feet at Mach 12. 
Time points are shown at four second intervals in the ascent path 
of figure 18. 

Ascent Path Optimization Without Sonic Boom Constraints 

The vehicle considered is an early, fully reusable shuttle con- 
figuration. The control parameterization and vehicle configura- 
tion are illustrated in figure 19. Initial trajectory optimization 
studies were performed to establish the maximum payload ascent 
path. This path is defined as: 

Max. (M^) 

subject to the terminal constraints 
= h^ - 303805 = 0 

t2 = Yf = 0 
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ALTITUDE, 1000 Meters (1000 Feet) 
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FIGURE 18 TYPICAL MAXIMUM GROUND OVERPRESSURE ON HIGH; 
DYNAMIC PRESSURE PATH. 
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FIGURE 19 COMPARISON OF VARIATION AND PARAMETERIZED 
CONTROL HISTORIES 
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24495 = 0 



where = final mass 

= final altitude 
= final flight path angle 
= final velocity 

The problem was solved using the two alternative optimization pro- 
cedures contained within the atmospheric flight path optimization 
program ATOP; these are: 

1. Variational calculus 

2. Multivariable search 

■V 

Results of the trajectory optimization are presented in figure 
20 in the Mach altitude plane. The region of overpressure con- 
straint violation is shown in the shaded area. The region is 
bounded on the left by the 2 psf overpressure boundary and on the 
right by the superboom (ray tube area goes to zero at the ground 
- theoretically infinite overpressure) . The constraint region 
is approximate and determined from a large number of sonic boom 
calculations at various flight conditions. The "optimum" trajec- 
tory passes through the region in the interval Mach 2 to Mach 3 
and 40000 to 70000 feet altitude. 

It is evident from figure 20 that the "optimxim" trajectory not 
only violated the sonic boom constraint but in accordance with 
the criteria presented, there is no feasible trajectory that will 
not violate it. However, the higher the altitude from which the 
boom emanates, the less severe the ground level superboom. This 
effect is not predicted by the method. Therefore, the trajectory 
optimization based on sonic boom constraint was conducted by rais- 
ing the superboom altitude as illustrated in figure 21. 

The resulting ascent paths are illustrated in the Mach altitude 
plane in figure 22. As the anticipated altitude at which the 
superboom boundary is crossed increased, orbital payload falls 
as indicated by the solid line in figure 23. This behavior was, 
therefore, consistent with the trajectories which were used to 
generate the approximate Mach altitude region in which ground 
overpressures in excess of 2 psf would be generated. 

A final optimal ascent was then computed in which the integral 
of the time in which a ground overpressure in excess of 2 psf was 
minimized while the orbital payload was constrained to three per 
cent less than the now known maximum value, that is. 
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ALTITUDE, 1000 FEET 


UPPER BOUNDARY 



FIGURE 20 APPROXIMATE OVERPRESSURE CONSTRAINT BOUNDARIES. 
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ALTITUDE, lOOOM (1000 FEET) 


IKUTfAL ASCENT FLIGHT PATH 



LOSS) 


FIGURE 22 


OPTIMAL TRAJECTORIES BASED ON SONIC BOOM CONSTRAINTS. 



TERMINAL MASS, KG (SLUGS) 



FIGURE 23 PERFORMANCE DEGREDATION AS A RESULT OF 
RAISING SUPERBOOM ALTITUDE. 
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Min (F 2 ) with Ap = 2 
f'l = hf - 303805 = 0 

'i'2 = Yf = 0 

- 24495 = 0 

- 10700 = 0 

Thus, in this case, the ground overpressure is minimized directly 
and no assumption is made regarding the location of the ground 
overpressure violation in the Mach altitude plane. 

The results are incorporated into figures 22 and 23. In figure 
22 the actual Mach altitude plane region responsible for ground 
overpressures in excess of 2 psf has moved slightly to the right 
of the assumed value. However, the duration of the violation is 
slightly less than anticipated as could be expected from the 
optimization process which properly accounts for all factors 
entering into the ground overpressure level not simply the Mach 
altitude effects. From figure 23 it can be seen that the cost 
in orbital payload required to raise the flight altitude at which 
a ground superboom is generated is approximately 1/2 per cent 
per 10,000 feet of altitude. 

A Note on the Ascent Solutions and the Superboom Problem 

Superbooms are created when the wave ray tubes emanated by the 
vehicle become very small. In the present report, the superboom 
boundary is taken to be the point at which a ray tube area goes 
to zero. This would indicate an infinite pressure in the over- 
pressure estimation methods now available. In actuality, the 
condition of a ray tube going to zero indicates a very high local- 
ized ground overpressure. In the studies of this section, this 
condition cannot be eliminated; however, the higher the vehicle 
altitude, which propagates a signature (which ultimately produces 
a ground superboom) , the less intense the resulting ground over- 
pressure will be. This was the rationale behind attempting to 
raise the anticipated superboom boundary crossing point. 

The question arises as to whether or not the superboom can be 
eliminated. This is equivalent to requiring one of two conditons 
to occur, either: 

1. The rays must still be coalescing at sea level. 
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2. The ray tube must go to zero at an altitude significantly 
greater than sea level. 

The space shuttle generates a continuous spectrum of signatures 
along its flight path, and barring severe nonstandard atmospheric 
discontinuities (or possible use of roll angle) , there is no 
manner in which the vehicle trajectory can suddenly pass from 
Type 1 to Type 2 or vice versa. 

Now, at higher altitudes of, say, 100,000 feet conditions leading 
to the ray tube going to zero at altitude are readily created. 

For example, in figure 20 the entire region of the supersonic Mach 
altitude plane to the left of the superboom boundary is such a 
region. As the trajectory approaches the superboom boundary from 
the left in this plane, the condition of a ray tube approaching 
zero at precisely ground level is approached. Once through 
the boundary, finite pressure signatures are experienced at the 
ground and the superboom is effectively underground. This behavior 
is illustrated in figure 24. 

In figure 24 (A) the vehicle lies to the left of the superboom 
boundary, and no sonic boom is felt at the ground. In 24 (B) the 
ground superboom occurs and the ray tube goes to zero at precisely 
sea level. In 24 (C) the vehicle lies to the right of the super- 
-boom boundary, and finite ground pressures are generated. 

Detailed examination of the factors entering into the generation 
of the superboom appears to indicate that a key factor in the ray 
tube focusing is the fact that the space shuttle trajectory is 
largely convex upwards. With this shape, after the first few 
moments of flight the flight path angle derivative, Y/ is negative. 
This trajectory behavior tends to lead to a focusing of the wave 
and the attendant superboom. 

Descent Paths 

Space shuttle orbiter aerodynamic characteristics and a nominal 
descent path pitch control history were supplied by NASA's Johnson 
Space Center. The nominal control history generated the reentry 
Mach altitude profile history of figure 25 in the ATOP program. 
Figure 26 presents the ground overpressure history produced by 
the nominal reentry. A maximum overpressure of 1.74 psf was 
encountered in this reentry path. The vehicle emitted the wave 
which propagated into this worst sea level condition while descend- 
ing through 65,000 feet at Mach 1.25. 

It can be seen from figure 26 that the space shuttle descent maxi- 
miam ground overpressure problem is essentially one of a temninal 
maneuver. Above a velocity of 4500 feet per second the maximum 
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OVERPRESSURE, Ap KPS (PSF) ALTITUDE, KM (FT) 



VELOCITY, KM/SECOND (FPS) 

FIGURE 25 NOMINAL REENTRY PATH 


9.76 



FIGURE 26 NOMINAL PATH OVERPRESSURES 
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overpressure is less than 1.0 psf. Accordingly, a one parameter 
family of descent paths characterized by constant angle of attack 
control was optimized to produce the minimum ground overpressure 
with trajectory cutoff at Mach 1.05. That is: 

Min (Fj^) 

Y. = (M-- - 1.05) = 0.0 

1 Nj 

The angle of attack range investigated was 0. < a < 60°^ Mini- 
miim boom was found at the maximvim angle of attack of 60 where 
the maximum ground overpressure was found to be 1.32 psf pro- 
pagated from 99,000 feet at Mach 1.88. This minimum "maximiam 
ground overpressure" path is presented in figure 27 together 
with the JSC nominal descent and two intermediate paths produc- 
ing ground overpressures of 1.62 psf at a = 24^ and 1.58 psf at 
a = 40°. Introducing the high angle of attack as the space 
shuttle vehicle passes through Mach 5 produces a high lofting 
maneuver which reduces the Mach number at a given altitude. The 
flight path point responsible for the maximum overpressure then 
tends to move to higher altitude and Mach number with increasing 
angle of attack while the magnitude of the maximum overpressure 
diminishes. 

It should be noted that the aerodynamic information provided by 
JSC is limited to a max = 24°. Thus, the program used extra- 
polated aerodynamics on the higher angle of attack reentry paths. 
The variation of maximum ground overpressure with the angle of 
attack value is presented in figure 28. If angle of attack is 
limited to 24°, then the corresponding maximum ground overpressure 
encountered is 1.62 psf, a seven per cent (7%) decrease in the 
nominal trajectory. 

Further optimization studies were then undertaken using two and 
four parameters to represent the control history as discussed in 
reference 8. The results are summarized in figure 29. If the 
angle of attack range investigated is limited to 0 < a < 24°, the 
best two parameter solution produces a maximum ground overpressure 
of 1.43 psf. The best four parameter solution produces 1.42 psf. 

If the angle of attack range investigated is increased to 0 < a < 
60°, the best two parameter solution is identical to the constant 
a = 60° solution of 1.32 psf. By using four parameters a slight 
improvement is produced, and the maximum ground overpressure 
recorded is 1.31 psf. 

Subsequently, the reentry paths were optimized starting from the 
point where the vehicle passed through Mach 10.0 using the angle 
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VELOCITY, KM/SECOND (EPS) 


FIGURE 27 


TERMINAL DESCENT PATHS, CONSTANT a AND NOMINAL. 













MAXIMUM OVERPRESSURE, Ap, KSM (PSF) 


9.76 



FIGURE 28 MAXIMUM OVERPRESSURES, CONSTANT « DESCENTS. 
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319NV M01 


319NV H9IH 


FIGURE 29 MAXIMUM OVERPRESSURES ON CONSTANT a REENTRIES, 





















of attack range 0 < a < 24°. The results are sioitimarized below. 
Only a slight gain resulted in going to angle of attack modula- 


tion at Mach 10 rather than 

Mach 5. 


TRAJECTORY 

APmax 

ANGLE OF ATTACK 

Best constant a 

1.53 

23.4° 

Best two parameter 

1.49 

9.7°, 23.9° 

Best four parameter 

1.45 

8°, 17°, 12.9°, 22°. 

Finally, some pitch control 
with the following results: 

optimization studies were undertaken 

TRAJECTORY 

APmax 

PITCH ANGLE 

Best constant pitch 

1.41 

12.04° 

Best two parameter 

1.40 

18.8°, 12.5° 


The last result was the best solution obtained without extrapolation 
of the MSC aerodynamic data. 
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USE OF THE PROGRAM 


Input to the program consists primarily of vehicular data, 
atmospheric properties, aircraft flight conditions and 
near-field pressure signatures data. Atmospheric properties 
are input in the form of temperature and wind profiles. 

The atmospheric pressure variations with altitude are com- 
puted in the program from input temperature profiles using 
the perfect gas law and the hydrostatic equation. Winds 
in both the northerly and easterly direction can be input 
allowing winds shear. However, vertical winds and atmospheric 
turbulence are not accounted for in this program. 

Flight conditions include Mach number, altitude, flight path 
angle and aircraft acceleration. The three components of 
aircraft acceleration are expressed in terms of the time 
rates of change of Mach number, flight path angle and head- 
ing respectively. The near-field pressure signature data 
consists of the signature itself and the corresponding 
location relative to the aircraft. The signature is input 
in the form of P/P (more commonly known as DP/P) versus 
X, where X is the special coordinate measured parallel to 
the aircraft velocity vector. The signature can be input 
directly and if done so, it should be consistent with the 
flight conditions specified. The flight conditions that 
affect the near-field signature are Mach nimber and lift 
coefficient. Therefore, the lift coefficient should first 
be estimated from the aircraft weight, flight altitude, 

Mach number and aircraft acceleration. The near-field 
signature corresponding to the Mach number and the required 
lift coefficient can then be determined experimentally by 
wind tunnel tests or by theoretical means. Normally a 
near-field signatures are determined in the wind tunnel. 
Signatures are obtained at several lift coefficients or 
angles of attack and in several azimuthal planes of the 
model. The near-field signature corresponding to the 
specific lift coefficient and locations relative to the 
model can then be estimated by interpolation. 

General Input Procedure 

The use of the SBOOM program is illustrated in figure 30. 
Measured pressure signatures from wind tunnel tests (or any 
other source) are stored in a pressure signature data base 
by the GETTAB program (see Volume II) . Vehicle properties 
and atmospheric conditions are read into the SBOOM program 
through the normal input channels. Flight conditions may 
also be read into the program by the same channels. However, 
flight conditions may be obtained from a trajectory program 
via an auxiliary input device (tape or disk). In figure 30, 
the POST program is illustrated. 
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MEASURED PRESSURE INPUT TRAJECTORY 

SIGNATURES FROM TflE CONTROL AND 
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FIGURE 30 ILLUSTRATION OF THE USE OF SBOOM. 








NAMELIST Data Format 


The program uses NAMELIST input for the following reasons; 

1. It is a simple name oriented input easily under- 
stood by most engineers. 

2. The format is standard and does not require 
relearning from program to program. 

3. It is easily modified by the engineer or programmer 
when adding input variables to the program. 

When NAMELIST read is encountered in a program, the entire 
input file is scanned up to an end-of-file or a record with 
a $ in column 2 followed immediately by the namelist name 
requested by the program. Succeeding data items are read 
until a second $ is encountered signifying the end of the 
NAMELIST. Any data on the input file before the requested 
namelist is found will be ignored. All data between the 
opening and closing $ are interpreted by NAMELIST. 

The data item within the NAMELIST statement may be in any 
of three forms; 

V = c, 

a = dj^ / . . . , d j , 
a(n) ' 

V is a variable name; c is a constant; a is an array name, 
and n is an integer constant subscript, d^ are simple con- 
stants or repeated constants of the form R*c, where k is 
the repetition factor. Data items and constants must be 
separated by commas. 

The number of constants, including repetitions, given for 
an unsubscripted array name must equal the number of elements 
in that array. For a subscripted array name, the number 
of constants need not equal, but may not exceed, the mamber 
of array elements needed to fill the array. 

The specified constant of the NAMELIST statement may be 
integer, real, double precision, complex of the form (c, ,Cp) 
or logical of the form T, or .TRUE., F, or .FALSE.. A 
logical or complex variable may be set only to a logical 
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and complex constant, respectively. Any other variable may 
be set to an integer, real or double precision constant. 

Such a constant is converted to the type of its associated 
variable. 

Constants and repeated constant fields may not include 
embedded blanks. Blanks, however, may appear elsewhere 
in data records. 

The entire card record excluding the first character is 
permitted. More than one card may be used for input data, 
and arrays may be split between cards. All except the last 
record must end with a constant followed by a comma, and 
no sequence numbers may appear. The first column of each 
record is ignored. 

The set of data items may consist of any subset of the 
variable names associated with the NAMELIST name. These 
names need not be in any particular order. 

SBOOM Program Input 

The input procedure for the SBOOM program is illustrated in 
figure 31. The initial read of the $SBIN namelist input 
establishes the vehicle properties and atmospheric data. 

The $SBIN input list is summarized in figure 32. The source 
for the flight conditions data is also established from the 
above read. After the flight conditions are determined, the 
corresponding pressure signature is determined by the geometry 
similarity rule. Then the Thomas extrapolation procedure is 
executed. The program then returns for new flight conditions. 

The flight conditions data can be obtained from four different 
sources depending on the input value of IROPTN read in the 
$SBIN namelist above. 

IROPTN = 0 Flight conditions read in the namelist $FCON. 

IROPTN = 1 Conditions read from a binary tape in fixed 

format. 

IROPTN = 2 Flight conditions read from a binary tape 
in the old POST format. 

IROPTN = 3 Flight conditions read from a binary tape 
in the new format (NAS 1-12165) . 
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READ $SBIN 
TAPE 5 


READ $FCON 
TAPE 5 


READ POST DATA 
TAPE 8* 


READ BINARY 
FORMATTED FILE 
TAPE 8 


INTERPOLATE/EXTRAPOLATE 
FOR PRESSURE SIGNATURE 
BY GEOMETRIC SIMILARITY RULE, 


CALCULATE SONIC BOOM 
BY THOMAS METHOD 


IROPTN = 2 OR 3 DEPENDING ON WHICH POST FORMAT IS EMPLOYED. 


FIGURE 31 SBOOM INPUT FLOW LOGIC 






NAMELIST 

Name 

Nominal 

Values 


PRINTN 

.FALSE. 

Printout switch for Thomas 
boom program 
FALSE - no printout 
TRUE - complete normal printout 

AL 

VEHLEN 

100 

Vehicle length in feet 

REEL 

REFLEC 

1. 9 

Ground reflection factor 

OCODE 

2 

1 - requested signature in milliseconds 

2 - requested signature in feet 

NTEMP 

7 

Number of atmospheric temperature 
points 

NWINDE 

2 

Number of easterly wind points 

NWINDN 

2 

Number of northerly wind points 

TO (I) 

100*0. 

Temperatures 

ZT (I) 

100*0. 

Altitude of temperature points in 
FT*1 o3 

ZWE (I) 

100*0. 

Altitude of easterly wind points 
in FT*1Q3 

VOE(I) 

100*0. 

Easterly wind velocities 

ZWN(I) 

100*0. 

Altitude of northerly wind points 
in FT*10^ 

VON(I) 

100*0. 

Northerly wind velocities 

NALT 

1 

Number of altitudes at which requested 
signature is desired. 

ALT(I) 

50*0. 

Altitude at which requested signa- 
ture is desired (feet) 


FIGURE 32A NAMELIST $SBIN INPUT. 
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NAMELIST Nominal 


Name 

Values 


INDSBM 

1 

Sonic Boom frequency indicator. 

ACCOVP 

to 

• 

o 

Maximum acceptable overpressure 
(psf) 

TIBOOM 

60. 

Flight condition data on tape is 
not processed before this time 
(trajectory time units) . 

T2BOOM 

160. 

Flight conditions data on tape 
will be ignored beyond the time 
(trajectory time units) . 

HRACON 

0. 

Altitude below which the ray tube 
constraint will be calculated (ft) . 

PLOT 

.FALSE. 

Plot option for pressure signature 


interpolation. If .TRUE., an input 
file for the independent plot program 
is generated which generates data 
base source signatures and the inter- 
polated signatures for each flight 
condition. 

.FALSE. Print option for pressure signature 

interpolation. If .TRUE. , the source 
signatures and interpolated signatures 
are printed. 

.FALSE. Print options for POST trajectory 

tape. If .TRUE. , tape will be dumped. 

0 Nxamber of lateral sonic boom extra- 

polations to be computed to the right 
and to the left of the ground track. 
(NSWEEP = 2 means 2 extrapolations 
each side . 

1*0 Steps in cross range (statute miles) 

to be used for computing lateral sonic 
boom extrapolations (overpressure 
observations written Unit 12. 


FIGURE 32B NAMELIST $SBIN INPUT. 


PRNTRP 

■POSTPT 

NSWEEP 

DELXRG 
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Input Flight Conditions . - If the read option (IROPTN = 0) is 
selected, the user supplies successive $FCON input lists for the 
desired flight conditions which he wishes to evaluate. The $FCON 
input list is summarized in figure 33. Many of the flight condi- 
tions have more than one name corresponding to (1) the original 
Thomas program names, (2) the ATOP input names and (3) the POST 
input names. Multiple flight conditions input is illustrated in 
figure 34. 

Flight Conditions Tape . - If IROPTN is set to 1, flight conditions 
data must be read in a binary format of the following twelve 
items on one record in order. 


TIME Trajectory Time (Sec) 

HGC7F Vehicle altitude (Feet) 

GAM7D Flight path angle (Deg) 

GAM7D1 Time derivative of flight path angle (deg/sec) 
SIG7D Heading angle (Deg) 

PHILD Latitude (Deg) 

THL7D Longitude (Deg) 

AMACH Vehicle Mach number 

AMACHl Time derivative of Mach number 

BA77D Bank of roll angle (Deg) 

SIG7D1 Time derivative of heading angle (Deg/sec) 

ALPHD Angle of attack (Deg) 

The above list of names and descriptions corresponds to the names 
from the ATOP .program. However, this same format may be generated 
from any trajectory program and stored on a binary file on the 
above format for use by the SBOOM program. If the flight condi- 
tions tape is used, the trajectory time range to be evaluated may 
be specified by the $SBIN input, TIBOOM and T2BOOM. 

Flight Conditions from POST . - Special read options (IROPTN = 2 
or 3) have been incorporated in the SBOOM program which peimits 
the regular output tape from the POST programs to be interrogated 
for the required flight condition information. The only limita- 
tion is that the first item on the POST output list for each print 
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NAMELIST 

Name 

MACH 

AMACH 

FLTALT 
HGC7F 
ALT I TO 

MOOT 
AMACH 1 
MACHDT 

PSIDOT 

SIG7D1 

GAMDOT 

GAM7D1 

FPA 

GAM7D 

GAMMAR 

LONGP 

THL7D 

LONG 

LATP 

PHILD 

GDLAT 

HEAD 

SIG.7D 

SIGMAR 

POG 

BA77D 

BNKANG 


ALPHA 

ALPHD 

TIME 


FIGURE 33A 


Nominal 

Values 

0. Vehicle Mach number 


100000. Vehicle Altitude (feet) 


0. 

Time derivative of Mach number 
(per sec) 

0. 

Time derivative of heading angle 
(deg/sec) 

0. 

Time derivative of flight path 
angle (deg/sec) 

0. 

Flight path angle (deg) 

0. 

Longitude (deg) 

0. 

Latitude (deg) 

0. 

Heading angle (deg) 

2116,2 

Atmospheric ground pressure (psf) 

0. 

Roll angle (deg) 

0. 

Angle of attack (deg) 

0. 

Trajectory Time (sec) 


NAMELIST $FC0N INPUT. 
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NAMELIST 

NAME 

NX 

X(I) 

DPP (I) 

ML 

ROVERL 

NALT 

ALT 

PHI 


FIGURE 33B 


NOMINAL 

VALUES 

0 

50*0 


50*0 


0 . 


0 . 


1 

50*0 

0 


Number of input values for X(I) and 
DPP (I) . 

Spatial coordinate measured parallel 
to the vehicle velocity vector, the 
independent variable for the near- 
field pressure signature (can be input 
in any unit of length) . 

The value of overpressure Ap/p correspond- 
ing to X (I) . 

Model reference length in the same units 
of length as for X(I). 

Radial distance (non-dimensional by the 
aircraft length) from the flight path 
corresponding to the input pressure 
signature . 

Number of altitudes at which extrapolated 
signatures are desired. 

The altitudes at which extrapolated 
signatures are desired. 

Angle of extrapolation measured from the 
vertical (downward) . The angle (PHI - 
BNKANG) must be the orientation of the 
signature used in the sonic boom calcula- 
tion. 


NAMELIST $FC0N INPUT. 
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$SBIN 


IROPTN = 0. 

I 

$ 



VEHICLE AND ATMOSPHERIC 
INPUT INCLUDING FLIGHT 
CONDITION OPTION 


FIRST FLIGHT CONDITION 


SECOND FLIGHT CONDITION 


PROGRAM TERMINATES ON 
END-OF-FILE 


FIGURE MULTIPLE FLIGHT CONDITIONS USING $FCON 



step must be TIME. The program interrogates the file for the 
other flight conditions specified in figure 35 until a new value 
of TIME is encountered. If all conditions are not present on 
the file, the nominal values shown in figure 35 are assxmied. 

It should be noted the POST output tapes are entirely different 
formats for IROPTN = 2 and IROPTN = 3. The user must be aware 
of which version of POST is being used. Option 3 uses the POST 
version generated under Contract NAS 1-12165. 

Basic Deck Setup . - The program may be used independently) 
sequentially with POST or other programs, or it may be used in 
ODIN for optimization or trajectory matching studies. The deck 
setup for the SBOOM program independent of a ancillary programs 
(see Volume II) consists of the above described program input data 
preceded by a set of control cards. The control cards vary depend- 
ing primarily upon the source of the flight conditions. The basic 
deck setup is illustrated in figure 36. The illustrated setup 
assumes the pressure signature data base and flight conditions 
have been previously stored on data cell. The file replacement 
parameters for the pressure signature data base and the flight 
conditions file are the fourth and fifth file parameter on the 
SBOOM execution card respectively. If the flight conditions are 
read from input or were generated by a prior execution in the 
same job, the FETCH card for FLTCON would not be necessary. The 
program loads in 61000 octal and executes in 53500 octal core 
locations. 

The file parameters for the SBOOM program include the following: 

NMLIST BCD output file for ODIN data base modification. 

FLTCON Binary flight conditions file. 

PSIGS Binary pressure signature data base file. 

SIGPLT BCD file for generating pressure signature inter- 

polations from the independent plot programs. 

FCONPLT Binary file of input flight conditions (from 

TIBOOM to T2BOOM) . 

CONTOUR Binary file of pressure signatures observations 
written when NSWEEP greater than zero. 

Deck Setup for Use of SBOOM with POST . - Figure 37 illustrates 
the use of SBOOM with POST for evaluation of sonic boom constraints 
along a POST generated trajectory. 

Deck Setup for Use of SBOOM within ODIN . - Figure 38 illustrates 
the use of SBOOM with ODIN for a sonic boom footprint study. 
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POST 

NAME 

NOMINAL 

■VAUJE.. 

DESCRIPTION 

TIME 

— 

TIME (sec) 

ALTITO 

0. 

ALTITUDE (ft) 

GAMMAR 

0. 

FLIGHT PATH ANGLE (deg) 

GAMDOT 

0. 

TIME DERIVATIVE OF FLIGHT PATH 
ANGLE (deg/sec) 

SIGMAR 

0. 

HEADING (deg) 

GDLAT 

0. 

LATITUDE (deg) 

LONG 

0. 

LONGITUDE (deg) 

MACH 

1.0101 

MACH NUMBER 

MACHDT 

0. 

TIME RATE OF MACH CHANGE (per sec) 

BNKANG 

0. 

BANK ANGLE (deg) 

SIGDOT 

0. 

TIME RATE OF CHANGE OF HEADING 
(deg/sec) 

ALPHA 

0. 

ANGLE OF ATTACK (deg) 


FIGURE 35 FLIGHT CONDITIONS FROM POST. 



JOB 

USER, 

FETCH , A3682 , SPRZ14 , B I NA RY , , FETCH . 
FETCH,DA423,,DATA,,,PSIGS. 

FETCH, A4190,, BINARY,, FLTCON. 
FETCH,A4193,, BINARY,, SBOOM. 
SBOOM,,,NMLIST,FLTCON,PSIGS,SIGPLT. 
7-8-9 

DATA 

6 - 7 - 8- 9 


FIGURE 36 DECK SETUP FOR SBOOFl USING STORED FLIGHT 
CONDITIONS TAPE. 
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JOB - - 


USER. 

FETCH , A3682 , SPRZ14 .BINARY,, FETCH . 
FETCH, DAA23, .DATA, ,,PSIGS. 

FETCH, A4199,,B I NARY,, POSTF. 

FETCH, A4193,, BINARY, .SBOOM. 

POSTF. 

SBOOM,,,,PROFIL. 

7-8-9 

(post data) 

7-8-9 

(SBOOM data) 

6-7-8-9 


FIGURE 37 DECK SETUP FOR USING POST WITH SBOOM. 



JOB - - - 
USER - - - 

FETCH , A3 6 8 2 , SPRZ 1 4 , B INARY , , GOGET . 

7-8-9 

7-8-9 

'CREATE DBASE' 

(INITIAL DATA BASE DEFINITIONS) 

*EOF 

'EXECUTE POST' 

(POST DATA) 

*EOF 

'EXECUTE VARIAN' 

(NO DATA - THIS EXECUTION INITIALIZES THE ODIN PLOT SYSTEM) 
*EOF 

'EXECUTE SBOOM' 

(SBOOM DATA) 

*EOF 

'EXECUTE PLOTTER' 

(PLOTTER DATA) 

*EOF 

'END ODIN' 


FIGURE 38 USE OF SBOOM WITHIN ODIN. 
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The loop is established around POST, SBOOM and PLOTTER whereby the 
flight conditions and trajectory performance are generated by POST, 
the sonic boom criteria are generated by SBOOM and the plotted 
results are generated by PLOTTER. The independent plot program 
is executed to generate the pressure signature contours along the 
flight path. (See reference 2 for complete details on the use of 
ODIN.) Appendix A discusses the results for the test case illustra 
ted in figure 38. 
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APPENDIX A - TEST CASE 

This appendix illustrates the type of analysis available from the 
expanded sonic boom prediction program by using it with the ODIN 
system and associated program modules. The analysis involved the 
use of the POST Program, the Sonic Boom Program and an independent 
plot program. POST was used for calculating and optimizing a 
reentry trajectory. Optimization was based on maximum cross 
range as the performance criteria and bank angle as a control 
function. Angle of attack was held constant through most of the 
reentry trajectory but was decreased in the linear manner during 
the latter portion of the reentry. The Mach-Altitude (M-h) 
flight profile for the trajectory is illustrated in figure A-1. 

The region in which the ground overpressure becomes a factor in 
the analysis is illustrated in the lower portion of the M-h 
plane below about Mach 3. Figure A-2 shows a time history of the 
control angles (angle of attack and bank angle) throughout the 
reentry. Angle of attack was fixed and bank angle was determined 
in the optimization process described above. 

Figure A-3 shows the ground track in terms of latitude and longi- 
tude. The region of ground overpressure significance is 
illustrated in the inset view, a small scale portion of the 
ground track near the end of the simulated trajectory. The out- 
line of the overpressure boundary (Ap/p > 0) on the ground track 
is illustrated in dash lines in the inset. Figure A-4 shows the 
sonic boom overpressure contours in terms of cross range and tra- 
jectory time. The zero overpressure boundaries of figure A-4 map 
into the overpressure boundary illustrated (by dash lines) in 
figure A-3. The contours illustrated in figure A-4 were generated 
by the independent plot program using a mesh of points and 
associated overpressure generated by the SBOOM Program. The 
sonic boom program has the capability of determining overpressure 
at various points on the ground by sweeping through a sequence of 
projection angles adjacent to the vertical plane. For each sweep 
angle the SBOOM Program extrapolates a sonic boom signature 
obtained from the data base onto the ground. Each sweep angle 
represents a point in the overpressure contour mesh. The result- 
ing mesh points from the extrapolations are stored for later use 
by the independent plot program in generating the contour map. 

It should be noted that the contour map does not represent a true 
projection of the sonic boom footprints, since it is plotted as 
a function of time along a curved ground track (and not in the 
latitude-longitude plane) . The peak overpressures near the lower 
portion of the figure represent focusing which is also indicated 
in the upper portion of the inset in figure A-3 . The data from 
figure A-4 has not been entirely mapped into the ground track of 
figure A-3. Cross cuts of maximum overpressure versus range from 
the ground track at fixed time are illustrated in figure A-5. 
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FIGURE A1 MACH-ALTITUDE PROFILE 




TIME (SECONDS) 

FIGURE A2 REENTRY CONTROL ANGLE HISTORIES 

















Peak overpressures are illustrated to the right of the ground 
track and represent the focusing illustrated in figures A-3 and 
A-4. 


The data presented herein represents an illustration of analyses 
which may be performed with the SBOOM Program when used within 
the ODIN system. Indeed the analysis can be extended using com- 
puter tools available within the ODIN system for studying the 
sonic boom problem. It should be noted that the present study 
does not represent a complete analysis.. For example, the focus- 
ing which produces high overpressures on the ground could be 
alleviated by "unbanking" the aircraft in the region where focus- 
ing occurs. Further, the data suggests that promising results 
could be obtained by “some trajectory optimizing on the basis of 
sonic boom overpressure rather than just cross range. 
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If Undeliverable (Section 158 
Po.stal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall he 
conducted so as to contribute . . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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technical information considered important, 
complete, and a lasting contribution to existing 
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TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 






